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Competing Spin Liquid States in the Spin- 1/2 Heisenberg Model On Triangular Lattice 


Wen-Jun Hu, Shou-Shu Gong,* Wei Zhu, and D. N. Sheng 
Department of Physics and Astronomy, California State University, Northridge, California 91330, USA 

We study the spin-1/2 Heisenberg model on the triangular lattice with the antiferromagnetic first (Ji) and 
second (J 2 ) nearest-neighbor interactions using density matrix renormalization group. By studying the spin 
correlation function, we find a 120° magnetic order phase for J 2 < 0.07Ji and a stripe antiferromagnetic 
phase for J 2 > 0.15Ji. Between these two phases, we identify a spin liquid region characterized by the 
exponential decaying spin and dimer correlations, as well as the large spin singlet and triplet excitation gaps on 
finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which 
indicates more than one spin liquid candidates in this region. While the sector with spinon is found to respect 
the time reversal symmetry, the even sector without a spinon breaks such a symmetry for finite-size systems. 
Furthermore, we detect the signature of the fractionalization by following the evolution of different ground 
states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we 
explore the nature of the spin liquid phase and find the optimal parameter region for the gapped Z 2 spin liquid. 

PACS numbers: 73.43.Nq, 75.10.Jm, 75.10.Kt 


Quantum spin liquids (SLs) are long-range entangled states 
with remarkable properties of fundamental importance [1]. 
The SL physics has been considered to be essential to under¬ 
stand strongly correlated systems and unconventional super¬ 
conductivity [2, 3]. The simplest and perhaps most striking 
SLs are the gapped topological SLs, which develop a topo¬ 
logical order [4-6] with the emergent fractionalized quasipar¬ 
ticles [7-9]. Although SLs have been studied intensively for 
two decades and demonstrated in contrived models [10-20], 
the microscopic condition for the emergence of SLs in frus¬ 
trated magnetic systems is not well understood. 

At the experimental side, possible SLs have been discov¬ 
ered in various materials. Among these materials, the most 
promising systems are the kagome antiferromangets includ¬ 
ing the Herbertsmithite and Kapellasite [21-25], as well as 
the organic Mott insulators with a triangular lattice structure 
such as k-(ET) 2 Cu 2 (CN )3 [26-29] and EtMe 3 Sb[Pd(dmit) 2]2 
[30, 31]. In all these materials, no magnetic order is observed 
at the temperature much lower than the interaction energy 
scale. These experimental findings have inspired intensive 
theoretical studies on the frustrated magnetic systems with 
strong frustration or competing interactions. 

Theoretically, the kagome Heisenberg model appears to 
possess a robust SL. Density matrix renormalization group 
(DMRG) studies suggest a gapped SL [32-35], which may 
be consistent with a Z 2 topological order [34, 35]. Varia¬ 
tional studies based on the projected fermionic parton wave 
functions however favor a gapless Dirac SL [36-38]. Interest¬ 
ingly, by introducing the second and third neighbor couplings 
[39-41] or the chiral interactions [42], DMRG [40-42] stud¬ 
ies recently discovered another topological SL — chiral spin 
liquid (CSL) [43, 44], which breaks time reversal symmetry 
(TRS) spontaneously and is identified as the ^ = 1/2 bosonic 
fractional quantum Hall state. On the other hand, the non¬ 
magnetic phases in the frustrated honeycomb and square Ji- 
J 2 models appear to be conventional valence-bond solid state 
[45-48]. 

The spin-1/2 triangular nearest-neighbor antiferromagnetic 


(AE) Heisenberg model was the first candidate proposed to re¬ 
alize a SL ground state by Anderson [2], although it turns out 
to still exhibit a 120° AE order [49-53]. To understand the tri¬ 
angular weak Mott insulator materials, combined theoretical 
and numerical studies [54-56] on a spin model with four-site 
ring-exchange couplings [57] find a gapless spin base metal 
with spinon Eermi sea. To enhance frustration [58-68], one 
way is to include the second-neighbor coupling J 2 , where a 
stripe ordered state is found with larger J 2 coupling [58, 59], 
and an intermediate non-magnetic region may emerge [60-63] 
based on different approaches. The variational Monte Carlo 
simulations find a nodal d-wave SL [61] and a gapless SL 
[62] as the candidates for this intermediate phase, although 
such method has the enhanced tendency of finding fraction¬ 
alized phase in nonmagnetic region. Very recently, a DMRG 
work [69] found the indication of a gapped SL which con¬ 
serves the TRS in the non-magnetic phase. However, the na¬ 
ture of the quantum phase with the intermediate J 2 remains 
far from clear. 

In this Letter, we study the spin-1/2 triangular model with 
the AE first and second nearest-neighbor Ji( Jj)-J 2 couplings 
based on DMRG calculations. The model Hamiltonian is 
given as 

i7 = Ji ^ S,-S,+J[ Y. S,-Sj+J2 Y 

(Fi)vertical (fij)zigzag ifd)) 

where the sums {i,j) and {{i,j)) run over all the first- and 
second-neighbor bonds, respectively. The first-neighbor cou¬ 
plings Ji and J[ are for the vertical and zigzag bonds as 
shown in Eig. 1(a). We study most systems with J( = Ji un¬ 
less we specify otherwise. We set Ji = 1 as the energy scale. 
By studying the spin correlations, we find a non-magnetic re¬ 
gion sandwiched by a 120° AE phase with three sublattices 
for J 2 < 0.07 and a stripe AE phase for J 2 > 0.15 as shown 
in Eig. 1. In this non-magnetic region, we identify two ground 
states with distinct properties in two sectors, indicating two 
competing candidates for SL phases. The spin and dimer cor¬ 
relations decay exponentially with small correlation lengths. 
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FIG. 1: (color online) Quantum phase diagram of the isotropic spin- 
1/2 J 1 -J 2 Heisenberg model on triangular lattice (Ji = J[). With 
growing J 2 , the system has a 120° AF phase for J 2 < 0.07, a stripe 
AF phase for J 2 > 0.15, and a SL phase in between. The schematic 
figures of the different phases also show the YC (a,b) and XC (c) 
cylinder geometries, (d)-(f) are the contour plots of spin structure 
factor for each phase. 


Interestingly, the chiral correlations decay exponentially fast 
for the ground state in the odd sector with an edge spinon, 
while it develops the long-range correlations in the even sec¬ 
tor (with no spinon) for finite-size systems, consistent with the 
level crossing between two SLs for the systems with different 
boundaries. The fractionalized spinon is detected through adi- 
abatically inserting spin flux. While the state in the odd sector 
agrees with a TRS preserving SL, the TRS breaking SL (e.g., 
chiral SL) may be a competing or nearby state in more ex¬ 
tended parameter space. Moreover, the strong anisotropy of 
bond energy along different directions is observed for some 
finite-size systems, which may imply a nematic order for 
gapped Z 2 SL [70]. This possible Z 2 SL is observed to be 
stabilized by a small bond coupling anisotropy (J[ > Ji), 
which suppresses chiral order in both sectors. 

We study the cylinder systems using highly accurate SU (2) 
DMRG [71, 72] for most of calculations and t/(l) DMRG 
[71] for inserting flux [40]. Two cylinder geometries known 
as XC and YC are studied, which have one lattice direction 
parallel to the x or y axis as shown in Fig. 1. We denote them 
as XCLy-Lx (YCLy-Lj;), where Ly and are the number of 
sites along the y and x directions, respectively. We study the 
cylinder systems with Ly up to 10 lattice spacings by keep¬ 
ing up to 20000 t7(l)-equivalent states in SU(2) DMRG and 
5000 states for inserting flux. The truncation errors are less 
than 10“® in all calculations, which leads to accurate results. 

Even and odd topological sectors .— Based on the resonat¬ 
ing valence-bond picture, the ground states of SL on cylinder 
can be either in the even or odd sector according to the parity 
of the number of bonds that are cut by a vertical line along the 



(b) J2=0.1 YC8-24, odd sector 



FIG. 2: (color online) The NN bond energy {Si ■ Sj) for J 2 =0.1 on 
the YC8-24 cylinder in (a) the even and (b) the odd sector. The left 
16 columns are shown here. The odd sector is obtained by removing 
one site in each boundary of cylinder. In both figures, all the bond 
energy have subtracted the average value —0.18. The red solid and 
blue dashed bonds denote the negative and positive bond energies 
after subtraction (with some numbers shown for clarity). 


enclosed direction. Usually, the odd sector can be obtained by 
removing or adding one site on each open edge of cylinder, 
which has been used successfully to find different sectors of 
the gapped SLs in kagome systems [33, 73, 74]. 

By doing simulations without and with removing one site 
at each boundary, we always find different ground states in 
these two sectors on YC cylinders {Ly = 6, 8,10), which are 
shown in Fig. 2 for YC8 cylinder as an example (See Sup¬ 
plemental Material for YClO cylinder) [75]. We find that the 
two sectors have the different bond energy distributions in the 
bulk of cylinder. While the vertical bonds are weaker in the 
even sector, they become the stronger ones in the odd sector 
as shown in Fig. 2. The nematic order, which is defined as 
the difference between the strong and weak bonds to describe 
the lattice rotational symmetry breaking, exhibits the distinct 
behaviors for the two states on our studied systems. While the 
nematic order grows with increased cylinder width in the odd 
sector, it decreases in the even sector [75]. 

Characteristic properties of different SL states .— Next, we 
further characterize the two states by studying correlation 
functions. In Fig. 3(a) we show the spin correlations for 
J 2 = 0.1, which decay faster with growing cylinder width in 
both states, indicating the vanishing spin order in both states. 
In Fig. 3(b), we demonstrate the dimer correlation function 
= {{Si ■ S,){Sk ■ Si)) - {Si ■ Sff{Sk ■ Si), which 
also decay exponentially to vanish. Interestingly, the states 
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FIG. 3: (color online) (a) and (b) are the spin and dimer correla¬ 
tion functions on different YC cylinders for J 2 =0.1. ^ denotes 
the corresponding correlation lengths on YClO cylinder in the odd 
sector, (c) Chiral correlations for J 2 = 0.1 on different cylinders. 
All the data are obtained from real number DMRG calculations ex¬ 
cept “YClO, even(c)”, which denotes the data obtained from com¬ 
plex wave function DMRG calculations on YClO cylinder in the 
even sector, (d) Chiral order from the boundary to the bulk for 
the bond anisotropic system in the even sector, which are obtained 
from the complex DMRG. (e) Schematic phase diagram for the bond 
anisotropic system at J 2 = 0.1. 


in the odd sector have very short correlation lengths almost 
independent of the system width. But in the even sector, the 
correlation lengths are longer than those in the odd sector for 
the smaller Ly = 6 and 8, which decrease with growing sys¬ 
tem width. 

To study the possible TRS breaking, we calculate the scalar 
chiral correlation function (xiXj) iXi = {Si^iX Si^ 2 )-In 
both sectors, the chiral order has the same pattern, where the 
up and down triangles have the same chirality direction. As 
shown in Fig. 3(c), in the odd sector, the chiral correlations 
decay quite fast to vanish without developing long-range or¬ 
der. However, the chiral order and spontaneous TRS breaking 
are very robust for the Ly = 6,8 systems in the even sector 
at the intermediate phase, where additional two fold ground 
state degeneracy is also obtained in the DMRG simulation for 
each system (with opposite chirality) associated with the TRS 
breaking. As we increase system width to Ly — 10, the chi¬ 
ral correlation becomes less robust, where different results are 
obtained depending on if we use complex or real initial wave 
function in DMRG simulation. The chiral correlation remains 
to be long-ranged in the complex wave function, where the 
TRS is spontaneously broken. However, if we use real num- 
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TABLE I: The bulk energy per site in tbe even (Ceven) and odd (eodd) 
sectors, the energy difference Ae = Ceven — Eodd, the spin triplet 
(At) and singlet (As) gaps in the odd sector for J 2 = 0.1 and 0.125 
on the YCLy (Ly = 6,8,10) cylinders. We use the fully converged 
results for Ly = 6 and 8, and the extrapolated results for Ly = 10 
as shown in Fig. 4. 


her wave function, the DMRG will find a state with short- 
range chiral correlations (the real state has near identical bulk 
energy as the complex state but higher energy near the edge). 

To further clarify the chiral order in the system, we consider 
a bond anisotropy perturbation by tuning the nearest-neighbor 
zigzag bond strength as J[ (see Fig. 1(a)). For J 2 = 0.1, we 
find that the SL region persists for 0.95 < J[ < 1.05. In the 
odd sector, the chiral order vanishes, and all the properties are 
consistent with J( = 1.0. In the even sector, the chiral order 
appears stronger for 0.95 < J'l % 0.99, where it grows a bit 
from YC8 to YClO cylinder (see Fig. 3(d) for J[ = 0.98). 
For J[ > 1.0, the chiral order decays pretty fast from the 
boundary to the bulk especially for the larger YClO cylinder, 
which indicates a possible vanishing of the chiral order in this 
region in the thermodynamic limit. At J[ = 1.0, the chiral 
correlations are strong and show long-range behavior, but the 
chiral order also decays with the increase of system width. 
Thus, the two states in both sectors may recover TRS at large 
system limit for 1.0 J'l < 1.05, which is the most possible 
region for stabilizing a Z 2 SL. On the other hand, for 0.95 < 
J'l 0.99, chiral order becomes stable in the even sector, 
while the fate of such a phase remains unclear depending on 
if the chiral order would develop in the spinon sector in the 
thermodynamic limit. We illustrate our finding in the phase 
diagram Fig. 3(e). 

We calculate the bulk ground-state energy in both sectors. 
The energy per site for different systems are presented in Ta¬ 
ble I for J 2 = 0.1,0.125. For the smaller system widths 
{Ly = 6 and 8), the odd sectors generally have the lower en¬ 
ergy than the even sectors, which lead to a positive energy 
splitting Ae = eeven — Codd- However, this splitting drops 
very fast with the increase of Ly, which is tiny for system 
width Ly = 10 (for example Ae ~ —0.0004 for J 2 = 0.1, 
see Fig. 4) indicating the close energy for states in both sec¬ 
tors. We also compute the singlet Ag and triplet A^ gaps 
by obtaining the ground state first and then sweeping the two 
low-lying states simultaneously or the 5 = 1 sector in the 
bulk of cylinder [67]. The results are shown in Table I. 

Inserting flux and the nature of different sectors .— Inserting 
flux is an effective way to find different sectors and to deter- 
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FIG. 4: (color online) Bulk energy per site e versus DMRG trunca¬ 
tion error e for the ground states in the even and odd sectors for (a) 
J 2 = 0.1 and (b) J 2 = 0.125 on the YC10-24 cylinder. The num¬ 
bers denote the kept 5(7(2) states. The energy data are fitted using 
the formula e(£) = e(0) -\- ae -\- he^. 


mine the basic properties of the quasi-particles in the ground 
states, which has been applied in DMRG to study different 
topological SLs [40, 73, 74]. To introduce a flux, we impose 
the twist boundary condition in the y direction by replacing 
terms Sf SJ + h.c. —)■ SJ + h.c. for all neighboring 
bonds crossing the y-boundary. We start from the even sector 
by adiabatically increasing 9 and measuring the evolution of 
the spin -2 local magnetization (S'f). With increasing 6, a net 
spin -2 accumulates on the open edges as shown in Fig. 5(a), 
which indicates that the quasiparticle responding to the flux 
here must carry spin, such as the spinon in chiral spin liquid 
[40] and the fermionic spinon (spinon bonded with vison) in 
Z 2 SL [73, 76, 77]. 

With the flux 0 = 0 —27r, the net spin grows continu¬ 
ously from 0 to 0.5 on the edges as shown in Fig. 5(b). At 
9 = 27r, an 5^ = ±1/2 spinon develops on each boundary, 
and the ground state evolves to a new sector. By further in¬ 
creasing 9 = 2Tr —i' Ttt, the net spin dissipates continuously 
to zero, and the system evolves back to the even sector. In 
Figs. 5(c) and 5(d), we demonstrate the entanglement spec¬ 
trum (ES) with inserting flux. At 6* = 0 in the even sector, 
the whole ES is symmetric about = 0. At 9 = 27r, the ES 
evolves to be symmetric about = 1 /2, which is consistent 
with the observed fractionalized spin-carrying quasiparticles 
on boundaries. At 0 = Ttt, the ES becomes the same as that 
at 9 = 0, indicating that the system evolves back to the even 
sector. Interestingly, each eigenvalue of the ES at 0 = 27r is 
doubly degenerate. By comparing the odd sector obtained by 
removing sites and the sector with inserted flux 9 = 2tt, we 
find that these two states have the same bulk energy and ES, 
indicating that the two sectors obtained by different methods 
are exactly the same. This odd sector might be consistent with 
the fermionic spinon sector of the Z 2 SL [76-78]. 

Summary and Discussions. — By means of DMRG calcula¬ 
tions on wide cylinder systems (up to Ly = 10 lattice sites) 
of the spin-1/2 J 1 -J 2 Heisenberg model on the triangular lat¬ 
tice, we find a SL region bordered by a 120° Neel AF phase 


(a) J 2 = 0 . 1 , YC 8 -I 6 , e=2n, AQr=0.48 



FIG. 5: (color online) (a) Real-space configuration of spin magneti¬ 
zation (5f) after adiabatically inserting a flux 6 = 27r. The red and 
blue circles denote the positive and negative (5f) with the area of 
circle proportional to the amplitude of (5f). AQr is the net spin ac¬ 
cumulation on the right edge AQr = where i denotes the 

site on the right edge of the cylinder, (b) Real-space configuration of 
the accumulated spin magnetization in each column with increasing 
flux. Low-lying ES at the flux (c) 9 — 0 and (d) 9 = 27r. The num¬ 
bers in ES denote the near-degenerate eigenvalues in the large weight 
levels. At 0 = 27r in (d), all the eigenvalues are double degenerate as 
explicitly shown by the circles for the low-lying levels. 


at small J 2 < 0.07 and a stripe AF phase for J 2 > 0.15. The 
spin and dimer correlations are shown to decay fast for wider 
systems with small correlation lengths comparable to lattice 
constant with large spin and singlet excitation gaps. The ES 
in the odd sector could be consistent with the theoretical de¬ 
scription of the fermionic spinon in the topological theory of 
Z 2 SL. However, the long-range chiral correlation is observed 
to be strong in the even sector for a space isotropic model. 
By tuning the anisotropic bond coupling, we find that the pos¬ 
sible gapped Z 2 SL is stabilized by some weak anisotropy 
(J( ~ 1.02). The chiral correlations are enhanced in the even 
sector by the opposite anisotropy (J[ ~ 0.98), which may be 
stabilized in both sectors by TRS breaking terms and we leave 
this open issue for future studies. 
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Note added .—While completing this work, we became 
aware of some related papers [69, 79]. We reached the similar 
conclusion on gapped SL with Ref. [69]. 
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Supplemental Material: Competing Spin Liquid States in the 
Spin-1/2 Heisenberg Model On Triangular Lattice 

I. Correlation functions 

First of all, we demonstrate the spin and dimer correlation 
functions for different J 2 on YC8-24 cylinder as shown in Fig. 
6 . In this figure, the correlation functions for J 2 = 0.1 and 
0.125 are shown in the odd sector. In Fig. 6 (a), the spin cor¬ 
relations I {Si ■ I decay slowly as a function of distance for 
J 2 = 0, 0.05, and 0.18, indicating the long-range magnetic 
order; instead for J 2 = 0.1 and 0.125, it decays exponentially 
to vanish, which is consistent with the absent magnetic order. 

In order to investigate the possible valence-bond solid or¬ 
der, we study the dimer-dimer correlation function on cylinder 
systems, which is defined as 

{^ij^kl) {^ij^{^kl)i ( 1 ) 

where {i,j) and {k,l) represent the nearest-neighbor (NN) 
bonds, and Bij = St ■ Sj. In Fig. 7 we show the dimer cor¬ 
relations for J 2 = 0.1 on YC8-24 cylinder in both the even 
and odd sectors. The dimer correlation in the odd sector de¬ 
cays faster than that in the even sector, which is shown in Fig. 
3(b) of the main text. In Fig. 6 (b), we show the dimer corre¬ 
lations for different J 2 , which all decay fast to vanish with an 
exponential manner. 

To detect the possible time-reversal symmetry (TRS) break¬ 
ing, we measure the chiral-chiral correlation functions for the 
four kinds of triangles Ai, A 2 , A 3 , and A 4 as shown in 
Fig. 8 (a) in both the even and odd sectors. The chiral-chiral 
correlation function is defined as 

{X^XJ) = ([^.,1 • (^..2 X S,,3)][Sj,i ■ {Sj ,2 X . 3 )]), (2) 

where Xi = Si^ ■ {Si ^2 x Si^s) is the scalar chiral order pa¬ 
rameter of triangle A^. As an example, we show the chiral 
correlations for J 2 = 0.1 on the YC8-24 cylinder in Fig. 8 . 
In the odd sector, all the chiral correlations decay quite fast 
to vanish (see Fig. 8 (c)), which is also observed on YC 6 and 
YClO cylinders and indicates no TRS breaking. However, in 
the even sector as shown in Fig. 8 (b), the chiral correlations 
for the triangles Ai, A 2 and A 3 exhibit the long-range chiral 
order, while those for A 4 decay fast. Using the complex num¬ 
ber DMRG, we can measure the chiral orders (xi) for all the 
triangles directly. In the odd sector, the chiral orders for all 
the triangles are vanished. In the even sector, we find that the 
chiral orders for Ai, A 2 and A 3 are finite, and (xi) for A 4 
is much smaller (for J 2 = 0.1 on YClO cylinder, it is about 
10“'*). In Fig. 8 (d) we show the chiral correlations of Ai by 
keeping different SU (2) DMRG states for J 2 = 0.1, J{ = 1.0 
and 0.98 on YClO cylinder in the even sector. By increasing 
the kept SU (2) DMRG states from 2000 to 4000, the chiral 
correlations change slightly, indicating the near convergence 
of the chiral correlations with increasing states. The non-zero 
chiral orders indicate the TRS breaking in the even sector. 



distance distance 


FIG. 6: (color online) (a) and (b) are the spin and dimer correlation 
functions for various J 2 couplings on the YC8-24 cylinder. In the 
spin liquid phase region, both correlation functions decay exponen¬ 
tially to vanish. 


II. Bond energy distributions in different sectors 

In Fig. 2 of the main text, we have shown the different NN 
bond energy distributions of the ground states in different sec¬ 
tors on YC 8 cylinder. Here, we demonstrate the bond energy 


(a) J2=0.1, YC8-24, even sector 




FIG. 7: (color online) The dimer-dimer correlation function for 
J 2 = 0.1 on YC8-24 cylinder in (a) the even sector and (b) the odd 
sector. The black bond in the middle is the reference bond. The blue 
and red bonds denote the positive and negative dimer correlations, 
respectively. 
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FIG. 8: (color online) (a) The four three-spin triangles which are 
used to calculate the chira-chiral correlation functions in (b) and (c). 
(b) and (c) are the chiral-chiral correlation functions versus triangle 
distance for J 2 = 0.1 on YC8-24 cylinder in the even and odd sec¬ 
tors. These chiral correlations are obtained from the real number 
DMRG. (d) is the chiral-chiral correlation functions of the triangles 
Ai at J 2 = 0.1 with anisotropic bond coupling J[ = 1.0 and 0.98 
in the even sector on YClO-20 cylinder with different kept SU{2) 
DMRG states. 


ing the same optimal states such as Msu{ 2 ) = 2000,3000, 
we find that the two states have the near identical bulk energy, 
but the real state has the higher total energy as there are much 
stronger bond energy fluctuations on the edge. This indicates 
that the real wave function is less converged near the bound¬ 
ary. 

In Fig. 10, we show the bond energy distributions obtained 
by keeping the most states for both the even and odd sectors, 
where the fluctuations are much smaller compared with the 
data in Fig. 9. The real and complex states in the even sector 
have the same pattern of the bond energy distribution in the 
bulk, which has the weaker vertical bonds. In contrast, the 
vertical bonds become the stronger bonds in the odd sector. 
This difference between the two sectors is consistent with our 
observations for YC8 cylinder, which are shown in Fig. 2 of 
the main text. 


III. Nematic order 

In Fig. 10, we notice that the NN bond energy distribu¬ 
tions have the lattice anisotropy with different bond energy 
along different lattice directions. In particular, the bond en¬ 
ergy in the odd sector seems to have a strong anisotropy as 
shown in Fig. 2(b) of the main text and Fig. 10(c) in Supple¬ 
mental Material. The strong lattice anisotropy may indicate 
a lattice rotational symmetry breaking. To investigate such a 
possibility, we study the anisotropy on different cylinder sys¬ 
tems. For convenience, we define the nematic order parame¬ 
ter (NOP) as the difference of the bond energy of the vertical 
bond and the others (NOP = Bond^igzag - Bondverticai)- In 
Fig. 11, we present the NOP for J 2 =0.1 and 0.125 on dif¬ 
ferent cylinders in two sectors, together with the results by 
tuning the anisotropic bond coupling J[. In the odd sector, 
the NOP increases continuously with growing system width 
for Ly = 6, 8,10 on YC cylinder for different J[. On the 
XC8 cylinder, the systems also have strong anisotropy. These 
results may suggest a finite NOP in thermodynamic limit. 
However, in the even sector the behaviors of the NOP ap¬ 
pear distinct for different J[. For = 1.0 and 1.02, the 
negative NOP on Ly = 6,8,10 indicates the weaker vertical 
bonds. However, for J[ = 0.98, the vertical bonds become 
the stronger bonds, and the NOP decreases continuously with 
growing system width for Ly = 6,8,10. 


distributions for J 2 = 0.1 on YClO cylinder in both sectors, 
and compare the results in the even sector obtained from the 
real and complex number states. As the even sector is harder 
to converge in DMRG calculations, we first show the bond 
energy distributions by keeping different states in the even 
sector. As shown in Fig. 9, by keeping 1000 SU{2) states, 
both the real and complex states exhibit some sort of string 
pattern. In particular, the real state shown in Fig. 9(a) has 
much stronger fluctuations. With increasing the kept states 
to 2000, the fluctuations become weaker. By comparing the 
ground-state energy of the real and complex states by keep- 


IV. Entanglement Entropy 

We have calculated the entanglement entropy (EE) on the 
Ly = 6, 8, 10, and 12 YC cylinders at J 2 = 0.1. In 
Eigs. 12(a)^(d), we show the entropy for Ly = 8, 10, and 12 
cylinders by keeping different 5(7(2) DMRG states Msjj( 2 ). 
As shown in Eig. 12(a)~(c), by using the less accurate data 
with the big truncation error, the linear extrapolated EE on 
YC8 cylinder and in the odd sector on YClO cylinder is far 
from the converged results. However, with the increasing kept 
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states, the trancation error decreases, and the more accurate 
EE is obtained. These three hgures suggest that it is hard to 
obtain the converged EE for the wider systems. In the even 
sector on YClO cylinder and odd sector on YC12 cylinder 
(Eig. 12(d)), we have kept up to 4000 SU{2) states, however, 
most data have the truncation error larger than 2x10 which 
indicates that they are far from convergency. By using these 
data, we only can obtain the lower bound of the EE with a lin¬ 


ear fitting, and due to the big truncation error, it is hard to get 
the reliable topological entanglement entropy on large cylin¬ 
ders. Based on the entropy data for Ly = 8 and 10 in the 
odd sector with the smaller error bar, we get the topological 
entanglement entropy close to In 2 from a reasonable fitting as 
shown in Eig. 12(e), which is consistent with the gapped Z 2 
spin liquid. 


11 


(a) J2=0.1 , YC10-24, Msu(2 )=1000, real, E=-122.390 



(c) J2=0.1, YClO-24, Msu(2 )=1000, complex, E=-122.425 
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(b) j2=0.1 , YC10-24, Msu(2 )=2000, real, E=-122.774 



(d) J2=0.1, YClO-24, Msu(2 )=2000, complex, E=-122.801 



FIG. 9: (color online) The NN bond textures for the even sector of J 2 =0.1 on the YClO-24 cylinder are calculated by using real [(a) and (b)] 
and complex [(c) and (d)] number DMRG. Different SU{2) states are kept: Msu{ 2 ) = 1000 in (a) and (c); Msu( 2 ) = 2000 in (b) and (d). 
The left 16 columns are shown here. In all figures, all the bond energy have subtracted a value —0.18. The red solid and blue dashed bonds 
denote the negative and positive bond textures, respectively. E is the total energy for each state. 
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(a) J2=0.1, YClO-24, Msu< 2 )= 5000 , real, E=-123.019 



FIG. 10: (color online) The NN bond textures for J 2 = 0.1 on YClO 
cylinders. The left 16 columns are shown here, (a) is the even sector 
obtained from real number DMRG by keeping 5000 5(7(2) states, 
(b) is the even sector obtained from complex number DMRG by 
keeping 3000 SU (2) states, (c) is the odd sector obtained by remov¬ 
ing one site in each boundary. In all figures, all the bond energy have 
subtracted a value —0.18. The numbers are the bond texture values. 
The red solid and blue dashed bonds denote the negative and positive 
bond textures, respectively. E is the total energy. For the system 
(a) by keeping 4000 5(7(2) states, the total energy E = —122.978, 
which is still higher than the energy of the complex state in (b). 
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FIG. 11; (color online) Cylinder width dependence of the nematic 
order parameter (NOP) for J 2 = 0.1 and 0.125 in (a) the odd sector 
and (b) the even sector. In the odd sector (a), the NOP increases 
with growing system width. In the even sector (b), we show the NOP 
obtained from both real and complex number DMRG calculations, 
which are consistent on YC6 and YC8 cylinders. On YClO cylinder, 
the complex states have the better convergence and thus exhibit the 
smaller fluctuations, (c) is the NOP in the even and odd sectors at 
J 2 =0.1 with tuning the bond anisotropy J[. 
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FIG. 12: (color online) (a)~(d) The entanglement entropy as a func¬ 
tion of the DMRG truncation error e at J 2 = 0.1 on YC8, YClO, and 
YC12 cylinders in the even and odd sectors. The largest kept SU (2) 
DMRG states Msu{ 2 ) are shown. In (a)~(c), we show the conver- 
gency of the EE by increasing the kept states. By doing the linear 
fitting with the less accurate data (red dashed lines), the extrapolated 
entropy is far from convergence. With more accurate data, we do the 
quadratic fitting (red solid lines), and the error bar is the difference 
between the extrapolated EE and the last raw data. In (d) we show 
the EE in the even sector on YClO cylinder and odd sector on YC12 
cylinder. Even with 4000 SU(2) DMRG states, most of these data 
have the truncation error larger than 2 x 10“®, which indicates that 
they are far from convergency. Based on the experience of (a)~(c), 
we do the linear fitting (dashed lines) to obtain the lower bound for 
the EE, and also make a quadratic extrapolation (solid lines) with the 
guessed big error bar for these two points, (e) The width dependence 
of the entanglement entropy in the even and odd sectors at J 2 = 0.1. 







